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Introduction 


Thermal  analysis  of  welds  seeks  to  determine  the  temperature  history,  i.e.,  change  in  temperature 
with  time,  at  positions  on  weld  cross  sections  [1,2].  Given  temperature  histories  as  a  function  of 
cross  section  position,  observed  weld  microstructure  and  weld  performance  can  be  correlated  with 
rates  of  temperature  change  during  welding,  and  residual  strains  within  welds  can  be  predicted  using 
models  for  which  temperature  histories  are  input  quantities.  The  direct-problem  or  first-principles 
approach  for  estimation  of  tempeature  histories  entails  solving  the  coupled  transport  equations  using 
numerical  methods  based  either  on  finite  differences  or  finite  volumes  [3,4],  This  approach  has  the 
advantage  of  providing  insight  concerning  the  relationship  between  weld  characteristics  and 
underlying  physical  processes,  i.e.,  the  nature  of  workpiece  and  energy  source  coupling.  For 
controlled  welding  experiments,  where  weld  meltpool  morphology  is  simple  and  thermocouple 
measurements  are  available,  the  direct  problem  approach  can  be  used  for  determination  of  material 
properties.  This  would  entail  adjustment  of  model  parameters  to  achieve  agreement  between 
predicted  and  measured  temperature  histories  at  thermocouple  locations.  The  direct  problem 
approach,  however,  has  many  disadvantages  for  practical  weld  thermal  analysis,  which  are  as 
follows. 

Direct-problem  based  models,  which  are  formulated  in  terms  of  basic  theory,  are  difficult  to 
apply  for  quantitative  weld  analysis  in  that  most  material  properties  needed  for  this  type  of  modeling 
are  not  available.  Weld  meltpool  morphologies  can  in  general  be  extremely  complex  due  to  the  wide 
variety  and  complexity  of  welding  processes,  e.g.,  hybrid  welding  processes.  Direct-problem  based 
models  are  not  structured  to  adopt  as  input  detailed  information  concerning  complex  shapes  of  weld 
meltpools.  That  is  to  say,  the  detailed  modeling  of  underlying  physical  processes  resulting  in 
complex  meltpool  morphologies  is  in  practice  not  feasible.  Direct-problem  based  models  are  not 
structured  to  adopt  thermocouple  measurements  as  input,  but  rather  only  as  verification  of  model 
output.  Finally,  direct-problem  based  models,  by  their  nature,  are  in  general  not  structured  for 
modeling  “overdetermined-systems”  [4],  i.e.,  weld  thermal  analysis  where  relatively  many 
experimental  measurements  are  available,  e.g.,  shape  features  of  solidification  boundaries  and 
thermocouple  measurements.  Accordingly,  with  respect  to  direct-problem  based  modeling, 
experimental  measurements  may  not  be  adaptable  as  model-input  information  or  be  useful  only  for 
model  verification. 

In  contrast  to  the  direct-problem  approach,  the  inverse-problem  approach  [5-7]  for  weld  thermal 
analysis,  using  parametric  functions  of  minimal  complexity,  has  many  advantages  for  quantitative 
estimation  of  weld  temperature  histories  [8-19].  This  approach  entails  calculation  of  temperature 
fields  within  the  volume  of  the  workpiece  using  parametric  functions  that  are  structured  for 
convenient  parameter  adjustment  with  respect  to  detailed  information  obtained  from  measurements, 
such  as  solidification  boundaries  and  thermocouple  measurements.  The  concept  of  system 
representation  using  “minimally  complex”  parametric  functions  is  related  to  that  of  system 
representation  using  linear  combination  of  basis  functions  [4],  which  follows  the  analysis  approach 
of  signal  processing,  where  analysis  is  in  terms  of  the  superposition  of  fundamental  modes  of  system 
response  [20],  Following  this  approach,  a  system’s  response,  no  matter  how  complex,  is  decomposed 
into  a  linear  combination  of  component  contributions  whose  formal  structure  are  characteristic 
modes  of  that  system.  For  example,  in  the  case  of  signal  processing  of  transmitted  waves,  e.g., 
electromagnetic  or  acoustic,  the  characteristic  modes  are  sine  and  cosine  functions,  which  represent 
in  principle  the  simplest  parameterization  in  terms  of  basis  functions  for  inverse  analysis.  Similarly, 
in  the  case  of  inverse  thermal  analysis,  the  characteristic  modes  are  the  kernel  and  fourier  series 
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solutions  to  the  heat  condunction  equation,  representing  the  simplest  parametrization  in  terms  of 
basis  functions  [21].  It  follows  that  the  inherent  complexity  of  any  numerical  procedure  for  inverse 
models  formulated  in  terms  of  linear  combinations  of  basis  functions  should  be  significantly  less 
than  those  for  applying  models  based  on  first  principles,  i.e.,  transport  equations,  such  as  finite- 
element  and  finite  volume  procedures  [3,4].  In  addition  to  being  physically  consistent  basis  functions 
for  parametric  representaion  of  weld  temperature  fields,  the  kernel  and  fourier  series  solutions  to  the 
heat  condunction  equation  have  general  functional  forms  that  are  sufficiently  flexible  so  as  not  to 
impose  any  form  of  model  bias  on  calculated  temperature  histories.  The  mathematical  foundation  of 
inverse  thermal  analysis  using  linear  combinations  of  basis  functions  is  that  of  least-squares 
parameter  optimization  [4,22,23], 

Inverse  analysis  compensates  for  lack  of  information  concerning  material  properties  and 
complex  physical  processes  associated  with  energy  source  and  workpiece  coupling  [24-27].  The  goal 
of  welding  models  formulated  in  terms  of  basic  theory,  inclusing  those  of  deep-penetration  welding, 
is  simulation  of  the  coupling  of  the  heat  source,  which  involves  melting,  fluid  flow  in  the  weld 
meltpool  and  heat  transfer  from  the  solidification  boundary  into  the  heat  affect  zone  (HAZ).  With 
respect  to  practical  analysis,  the  purpose  of  modeling  these  influences,  in  particular,  would  be 
generation  of  the  solidification  boundary,  the  surface  from  which  heat  is  transferred  into  the  HAZ, 
which  is  the  region  of  most  probable  weld  failure,  and  difficult  for  inspected.  Therefore,  inverse 
models  that  adopt  weld  cross  sections  of  solidification  boundaries  as  “model  input”  should  tend  to 
compensate  for  lack  of  information  concerning  physical  processes  occuring  within  the  melt  pool. 
This  reasoning  also  applies  to  physical  processes  occuring  within  the  mushy  zone.  With  respect  to 
inverse  analysis,  measured  solidification  boundaries  represent  estimates  of  mushy-zone  edges,  which 
are  at  the  solidus  temperature. 

Presented  here  is  a  case  study  inverse  thermal  analysis  of  structural  steel  deep-penetration 
welds.  This  analysis  provides  a  parameterization  of  temperature  histories  for  prediction  of  properties 
within  the  HAZ  of  welds  for  the  regime  considered.  The  present  study  applies  an  inverse  thermal 
analysis  procedure  that  uses  volumetric  constraints  on  calculated  three-dimensional  solidification 
boundaries.  These  constraints  are  such  that  two-dimensional  projections  of  calculated  solidification 
boundaries  are  constrained  to  fall  within  measured  solidification  cross-section  boundaries.  For  the 
present  study,  which  considers  deep-penetration  welds  of  structural  steel,  inverse  thermal  analysis 
includes  prediction  of  temperature  histories  at  experimentally  estimated  HAZ-edges.  The  parametric 
temperature  histories  presented  in  this  study  can  be  adopted  for  inverse  thermal  analysis  of  welds 
corresponding  to  other  process  parameters  or  welding  processes  whose  process  conditions  are  within 
similar  regimes.  Further,  this  study  can  contribute  to  a  parameter  space  for  inverse  themal  analysis, 
which  contains  a  wide  range  of  parameters  corresponding  to  different  welding  processes,  process 
conditions  and  different  types  of  metals  and  their  alloys. 

The  inverse  analysis  procedure  applied  here  entails  parametric  representation  of  the  temperature 
field  in  terms  of  numerical-analytical  basis  functions  [24-27],  which  consist  of  numerical  integrals  of 
analytic  functions  over  time.  The  conceptual  foundation  of  this  procedure  was  introduced  in 
reference  [28].  The  parametric  representation  provides  a  general  reduction  of  model  complexity  for 
purposes  of  weld  thermal  analysis.  Reduction  of  model  complexity  is  achieved  by  adopting 
numerical-analytical  models  for  either  the  heat  source  or  temperature  field  (or  both)  in  combination 
with  numerical  methods.  Among  studies  applying  this  approach  are  [2,  29-35], 

The  present  study  contributes  to  a  parameter  space  relating  shapes  of  weld  solidification 
boundaries,  or  shapes  of  any  specified  boundaries  within  the  weld,  to  weld  temperature  histories,  as 
well  as  to  weld  process  parameters.  The  inverse  analysis  methodology  applied  here,  which  permits 


2 


convenient  adjustment  of  parameters,  uses  discrete  source  distributions.  These  source  distributions, 
consisting  of  a  finite  number  of  point  sources  distributed  in  three-dimensions,  are  conveniently 
adjustable  with  respect  to  location  and  source  strength  for  calculation  of  temperature  fields  on  and 
exterior  to  specified  isothermal  boundaries,  e.g.,  solidification  boundaries  (see  reference  [24-27]). 
The  use  of  discrete  source  distrubutions  in  this  study,  which  are  within  three  dimensions  (i.e., 
volumetric),  is  formally  equivalent  to  the  Myhr-Grong  method  [30],  which  typically  employs 
discrete  source  distributions  within  two-dimensional  planes  transverse  to  the  relative  motion  of  heat 
source  and  workpiece.  References  [34,35],  and  reference  therein,  demonstrate  the  use  of  discrete 
source  distributions  for  numerical  modeling  of  transport  phenomena  associated  with  welding 
processes. 

An  important  part  of  the  present  study,  which  extends  the  methodology  for  inverse  thermal 
analysis,  is  examination  of  a  procedure  accessing  the  consistency  of  calculated  temperature  histories 
for  steel  welds,  which  uses  experimentally  measured  estimates  of  the  HAZ  edge.  In  order  to 
demonstrate  practical  application  of  this  procedure,  as  well  as  its  consistency,  a  reasonable  sampling 
of  well  characterized  steel  welds  (i.e.,  8  different  welds)  are  considered,  for  which  solidification  and 
estimated  HAZ-edge  boundaries  where  determined  experimentally  [36], 

The  organization  of  the  subject  areas  presented  are  as  follows.  First,  a  brief  description  of  the 
general  procedure  for  inverse  analysis  of  heat  deposition  processes  is  presented.  Second,  results  of 
inverse  thermal  analyses  of  structural  steel  deep  penetration  welds  are  presented.  These  results 
provide  a  quantitative  parametric  representation  of  temperature  histories  for  these  welds  and  for  any 
welds  associated  with  similar  welding  process  conditions.  Third,  a  discussion  is  presented  of  the 
results  and  of  using  experimentally  measured  HAZ  edges  for  a  consistency  check  of  calculated 
temperature  histories.  Finally,  a  conclusion  is  given. 

Inverse  Analysis  Procedure 

Following  the  inverse  analysis  approach,  a  parametric  model  provides  a  means  for  the  inclusion 
of  information  concerning  the  physical  characteristics  of  a  given  energy  deposition  process.  A 
physically  consistent  parametric  representation  of  temperature  fields  for  heat  deposition  during 
welding  of  plate  structures  is  given  by 


Nt  N, 

T(x,t )  =Ta  +^^C(xk)G(x,xk,K,nAt,V) 

k  — 1  n= 1 


and 


t(xio=t: 


where 


G(x, xk ,  t ,k, V)  =  -  exp 


{x-xk-vt?Hy-yk)2 

AKt 


r  oo 

2  2 

) 

1  +2^  exp 

[  m= 1 

KM  JT  t 

COS 

MJTZ 

COS 

MJTZk 

l 

i2 

i 

i 

] 

Nk 

and  C(x)=^Q(xk)d(x-xk). 

k= 1 


(Eql) 

(Eq2) 


(Eq  3) 


(Eq  4) 


3 


where  Q(xk )  is  the  value  of  the  discrete  source  function  at  location  xk.  The  quantities  k  ,  V  and  /  are 
the  thermal  diffusivity,  welding  speed  and  plate  thickness,  respectively.  The  constraint  conditions 
defined  by  Eq.(2),  which  represent  input  for  the  model  defined  by  Eqs.(l)-(4),  are  imposed  on  the 
temperature  field  by  minimization  of  the  objective  function  defined  by 

zr  =  (Eq5) 

where  Tnc  is  the  target  temperature  for  position  xcn  =  ( xcn,ycn,zcn ).  The  quantities  wn  ( n=l,...,N )  are 
weight  coefficients  that  specify  relative  levels  of  influence  associated  with  constraint  conditions  rnc . 
The  output  quantity  of  the  model  defined  by  Eqs.(l)-(4)  is  the  three-diemsional  temperature  field 
T(x,t )  spanning  the  entire  volume  of  the  workpiece. 

The  procedure  for  inverse  thermal  analysis  defined  by  Eqs.(l)-(5)  entails  adjustment  of  the 
parameters  C(xk)  ,  xk  and  At  .  The  parametric  model  combines  numerical  integration  with 
optimization  of  linear  combinations  of  numerical-analytical  basis  functions,  which  include 
fundamental  solutions  to  the  heat  conduction  equation  and  their  Fourier-series  representation  [21],  In 
particular,  Eq.(l)  defines  a  discrete  numerical  integration  over  time,  where  the  time  step  At  is 
specified  according  to  the  average  energy  deposited  during  the  time  At,  for  transition  of  the 
temperature  field  to  steady  state.  It  should  be  noted  that  the  formulation  of  the  inverse  analysis 
methodology  defined  by  Eq.(l)-(5)  is  equipped  with  a  mathematical  structure  that  satisfies  all 
boundary  conditions  associated  with  welding  of  plate  structures  (see  [24-27]  for  further  discussion). 

In  addition  to  the  parameters  defined  with  respect  to  Eqs.(l)-(5),  the  parameteric  model  applied 
for  inverse  analysis  includes  a  length  scale  parameter  Is,  where  in  general  Is  <  l  defined  by  Eq.(3),  for 
specification  of  the  spatial  scale  of  the  calculated  temperature  field  with  respect  to  which  parameters 
are  adjusted.  This  length  scale  parameter  provides  for  inclusion  of  more  details  of  shape  features  of 
measured  solidification  boundaries  to  be  adopted  as  constraint  conditions.  Inverse  analyses  using 
this  parameter  are  given  in  reference  [27,  37].  However,  for  the  analyses  that  follow,  the  two  length 
parameters,  i.e.,  the  plate  thickness  /  and  the  depth  Is  of  the  speciifed  region  of  the  temperature  field 
to  be  calculated,  are  taken  to  be  the  same. 

Table  1  Experimentally  measured  estimates  of  locations  (yc,zc)  on  solidification  and  HAZ-edge 
boundaries,  which  are  on  transverse  cross  section  of  Weld  1. 


SOLIDIFICATION 

ESTIMATED 

BOUNDARY 

HAZ-EDGE  i 

(yc  mm,  zc  mm) 

(yc  mm,  zc  mm)  , 

(5.052,  0.2105) 

(6.315,  0.2105) 

(3.579,  2.105) 

(5.684,  2.105) 

(2.316,  3.158) 

(5.052,3.158) 

(2.316,  4.210) 

(4.21,4.21) 

(1.895,  6.315) 

(3.579,  6.315) 

(1.895,  8.420) 

(3.368,  8.42) 

(1.684,  10.525) 

(3.158,  10.525) 

(2.316,  12.0) 

(3.789,12.0) 
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Table  2  Experimentally  measured  estimates  of  locations  (yc,zc)  on  solidification  and  HAZ-edge 
boundaries,  which  are  on  transverse  cross  section  of  Weld  2. 


SOLIDIFICATION 

ESTIMATED 

BOUNDARY 

HAZ-EDGE 

(yc  mm,  zc  mm) 

(yc  mm,  zc  mm)  \ 

(7.578,  0.2105) 

(9.262,0.2105) 

(5.894,  1.053) 

(8.42,  1.053) 

(4.421, 2.105) 

(7.578,2.105)  ! 

(3.158,  3.158) 

(6.736,  3.158) 

(2.316,  4.21) 

(5.473,4.21) 

(2.105,  5.263) 

(4.421,  5.263) 

(1.684,  7.368) 

(3.368,  7.368) 

(1.684,  9.473) 

(3.158,  9.473) 

(1.474,  10.525) 

(3.158,  10.525) 

(2.105,  12.0) 

(3.579,  12.0) 

Table  3  Experimentally  measured  estimates  of  locations  (yc,zc)  on  solidification  and  HAZ-edge 
boundaries,  which  are  on  transverse  cross  section  of  Weld  3. 


SOLIDIFICATION 

ESTIMATED 

BOUNDARY 

HAZ-EDGE 

(yc  mm,  zc  mm) 

(yc  mm,  zc  mm) 

(5.164,  0.2459) 

(7.377,  0.2459) 

(3.934,  2.459) 

(6.393, 2.459)  | 

(1.967,  4.918) 

(4.918,4.918) 

(2.213,  7.377) 

(3.443,  7.377) 

(1.721, 9.836) 

(3.197,9.836) 

(1.721,  12.30) 

(2.705,12.30)  I 

(1.475,  14.75) 

(2.705,  14.75) 

Table  4  Experimentally  measured  estimates  of  locations  (yc,zc)  of  solidification  and  HAZ-edge 
boundaries,  which  are  on  transverse  cross  section  of  Weld  4. 


SOLIDIFICATION 

ESTIMATED 

BOUNDARY 

HAZ-EDGE 

(yc  mm,  zc  mm) 

(yc  mm,  zc  mm) 

(5.090,  0.2679) 

(7.501,0.2679) 

(5.090,  1.340) 

(6.965,1.340) 

(4.019,  2.679) 

(6.430,  2.679) 

(2.947,  4.019) 

(5.626,4.019) 

(1.875,  5.358) 

(4.822,  5.358) 

(1.875,  8.037) 

(3.483,  8.037) 

(1.607,  10.72) 

(2.947,10.72) 

(1.875,  12.06) 

(2.947,12.06) 

(1.607,  13.40) 

(2.679,13.40) 

(1.607,  14.73) 

(2.679,14.73) 
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Table  5  Experimentally  measured  estimates  of  locations  (yc,zc)  on  solidification  and  HAZ-edge 
boundaries,  which  are  on  transverse  cross  section  of  Weld  5. 


SOLIDIFICATION 

ESTIMATED 

BOUNDARY 

HAZ-EDGE  ! 

(yc  mm,  zc  mm) 

(yc  mm,  zc  mm) 

(7.501, 0.2679) 

(10.18,  0.2679) 

(6.965,  1.340) 

(9.912,  1.340) 

(5.894,  2.679) 

(9.644,  2.679) 

(4.554,  4.019) 

(8.573,  4.019) 

(2.947,  5.358) 

(7.501, 5.358) 

(2.411, 6.698) 

(5.894,  6.698) 

(2.411, 8.037) 

(4.554,  8.037) 

(2.143,  9.377) 

(4.019,  9.377) 

(2.143,  10.72) 

(3.751,  10.72) 

(1.875,  13.40) 

(3.483,  13.40) 

(1.875,  14.73) 

(3.483,  14.73) 

Table  6  Experimentally  measured  estimates  of  locations  (yc,zc)  on  solidification  and  HAZ-edge 
boundaries,  which  are  on  transverse  cross  section  of  Weld  6. 


SOLIDIFICATION 

ESTIMATED 

BOUNDARY 

HAZ-EDGE 

(yc  mm,  zc  mm) 

(yc  mm,  zc  mm) 

(8.773,  0.2830) 

(10.47,  0.2830) 

(8.207,  1.415) 

(10.47,1.415) 

(6.792,  2.830) 

(9.905,  2.830) 

(5.377,  4.019) 

(9.056,4.019)  I 

(3.679,  5.66) 

(7.924,  5.66) 

(2.264,  7.075) 

(6.509,  7.075) 

(2.264,  8.037) 

(5.094,  8.037) 

(1.981,  11.32) 

(3.962,  1 1 .32) 

(1.981,  12.74) 

(3.962,  12.74) 

(2.264,  14.15) 

(3.962,  14.15) 

(2.264,  15.0) 

(3.962,  15.0) 

Table  7  Experimentally  measured  estimates  of  locations  (yc,zc)  on  solidification  and  HAZ-edge 
boundaries,  which  are  on  transverse  cross  section  of  Weld  7. 


SOLIDIFICATION 

ESTIMATED 

BOUNDARY 

HAZ-EDGE 

(yc  mm,  zc  mm) 

(yc  mm,  zc  mm) 

(9.793,  0.3159) 

(12.32,0.3159) 

(8.529,  1.580) 

(11.69,1.580) 

(6.634,  3.159) 

(10.11, 3.159) 

(4.739,  4.739) 

(8.845,  4.739) 

(3.159,  6.318) 

(6.95,6.318)  ! 

(2.211, 7.898) 

(5.686,  7.898) 

(2.843,  9.477) 

(4.423,  9.477) 

(2.843,  1 1 .06) 

(4.107,  11.06) 

(2.211,  12.64) 

(3.475,12.64)  j 

(1.589,  14.22) 

(3.159,  14.22) 

(2.211,  15.80) 

(3.475,  15.80) 

(2.527,  17.37) 

(3.791,17.37) 

(2.211,  18.0) 

(3.791,18.0)  | 

Table  8  Experimentally  measured  estimates  of  locations  (yc,zc)  on  solidification  and  HAZ-edge 
boundaries,  which  are  on  transverse  cross  section  of  Weld  8. 


SOLIDIFICATION 

ESTIMATED 

BOUNDARY 

HAZ-EDGE  ! 

(yc  mm,  zc  mm) 

(yc  mm,  zc  mm) 

(7.938,  0.3175) 

(12.38,  0.3175) 

(7.62,  3.175) 

(12.38,  3.175) 

(6.668,  4.763) 

(11.75,4.763)  ; 

(5.398,  6.35) 

(10.80,6.35) 

(4.128,  7.938) 

(9.842,  7.938) 

(3.81, 9.525) 

(8.255,  9.525) 

(2.858,  12.7) 

(6.35,12.7) 

(2.54,  15.89) 

(5.398,  15.89) 

(3.175,  19.05) 

(5.715,19.05) 

(2.858,  20.0) 

(6.033,  20.0) 
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Case  Study  Analysis  of  Steel  Welds 


In  this  section  results  of  inverse  thermal  analyses  of  structural  steel  deep-penetration  welds  are 
described,  which  correspond  to  different  weld  process  conditions  and  associated  process-control 
parameters.  The  significance  of  the  inverse-problem  approach  for  these  analyses  is  that  the  nature  of 
the  coupling  of  the  energy  source  to  the  workpiece,  which  is  a  function  of  beam  power  and  process 
control  parameters,  is  in  principle  difficult  to  specify  relative  to  analysis  based  on  the  direct-problem 
approach.  Previous  studies  considered  inverse  thermal  analysis  of  steel  welds  using  different  types  of 
constraint  conditions  [25,26],  The  present  study  uses  experimentally  estimated  solidification 
boundaries  for  assigning  volumetric  constraints  (see  Eq.(2))  on  the  calculated  temperature  fields.  As 
discussed  in  reference  [27],  for  regions  within  the  workpiece  not  close  to  the  energy  source, 
reasonable  estimates  of  temprature-field  values,  at  specified  locations,  are  sufficient  for  imposing 
constraint  conditions. 

For  steel  welds,  especially  those  of  structural  steel,  the  HAZ  consists  by  various  regions  whose 
microstructures  have  been  well  characterized  (see  reference  [1]).  These  regions  are  described 
schematically  in  Fig.  1.  Referring  to  this  figure,  it  should  be  noted  that,  for  inverse  thermal  analysis 
of  structural  steel  welds,  one  should  be  able  to  associate  the  observed  edge  of  the  HAZ  with  a  range 
of  temperatures  characteristic  of  that  region  of  the  weld  (see  Fig.  1). 


Fig.  1  Schematic  representation  of  different  regions  within  the  HAZ  for  steel  welds  and  approximate 
location  of  experimentally  observable  HAZ  edge. 
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Discrete  Source 


Top  Surface  Plane 


•  •  •  •  • 


Symmetry  Line 
at  Top  Surface 


Symmetry  Plane 


Fig.  2  Indexing  scheme  for  relative  locations  of  discrete  sources  C(xk),  k=l,...,Nk. 


The  structural  steel  deep-penetration  welds,  whose  inverse  analysis  is  presented  here,  consist  of 
laser  and  laser-GMA  hybrid  welds  [38].  The  analyses  presented  here  entail  calculation  of  the  steady 
state  temperature  field  for  a  specified  range  of  sizes  and  shapes  of  inner  surface  boundaries  defined 
by  the  solidification  boundary,  and  experimentally  observed  estimates  of  the  HAZ  edge.  The  shapes 
of  these  boundaries  are  determined  experimentally  by  analysis  of  transverse  weld  cross  sections 
showing  microstructure  revealing  solidification  and  estimated  HAZ-edge  boundaries.  For 
calculations  of  the  temperature  field,  which  adopt  solidification  boundaries  as  constraints,  the 
parameter  values  assumed  are  k=  5.88  x  10"6  mV 1 ,  Tm  =  1503.0  °C.  As  discussed  previously  [25, 
27],  reasonable  estimates  of  k  and  Tm  are  sufficient  for  inverse  analysis.  This  assumption  is 
sufficient,  within  reasonable  estimates,  in  that  the  set  of  parameters  C(xk),  k=l,...,Nk,  and /rare  not 
uniquely  determined  by  inverse  analysis.  Thus,  changing  estimated  values  of  k  would  require 
different  values  of  C( xk)  in  order  to  satisfy  specified  constraint  conditions  associated  with  Tm-  With 
respect  to  inverse  analysis,  the  interpretation  of  k  as  both  an  estimated  material  property  and 
adjustable  parameter  is  emphasized  within  the  following. 

The  goal  of  the  present  analysis  is  determination  of  a  set  of  parameters  that  can  serve  as  initial 
estimates  for  parameter  adjustment  with  respect  to  deep  penetration  welds  of  steels,  whose  process 
parameters  are  within  similar  regimes.  Parameter  adjustment  with  respect  to  other  welds,  which 
assume  the  results  of  this  study  as  initial  estimates,  would  adopt  k  and  Tm  as  adjustable  parameters, 
as  well  as  the  discrete  source  function  C(xk) .  Values  of  the  workpiece  thickness  /  and  welding  speed 
V  defined  in  Eq.  (3)  are  given  in  the  figures  below.  The  upstream  boundary  constraints  on  the 
temperature  field,  Tc  =  TM  for  (yc,zc)  defined  in  Eq.  (2),  are  given  in  Tables  1  through  8  for  the 
solidification  boundaries.  Given  in  Tables  9  through  16  are  values  of  the  discrete  source  function 
that  have  been  calculated  according  to  the  constraint  conditions  and  weld  specifications  given  in 
Tables  1  through  8.  The  relative  location  of  each  discrete  source  is  specified  according  to  indexing 
scheme  shown  in  Fig.  2.  Shown  in  Figs.  3  through  42  are  experimentally  measured  transverse  weld 
cross  sections  of  solidification  and  estimated  HAZ-edge  boundaries  [38],  and  different  planer  slices 
of  the  steady  state  temperature  field  that  have  been  calculated  according  to  the  constraint  conditions 
given  in  Tables  1  through  8  for  the  solidification  boundary.  Referring  to  the  planar  slices  of  the 
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calculated  temperature  fields  shown  in  these  figures,  it  can  be  seen  that  all  boundary  conditions  are 
satisfied,  namely  the  condition  T(x,t)-TM  at  the  solidification  boundary,  and  V7  • «  =  0  at  surface 

boundaries,  where  h  is  normal  to  the  surface. 

Table  9  Source  function  C(xk)  calculated  according  to  solidification-boundary  constraint  conditions 
given  in  Table  1,  where  A l  =  (12/60)  mm,  Xk  =  yk=  0.0  for  k  =  1  to  13  (WELD  1). 


k 

Q/20.0 

Zk  (A  / ) 

1 

0.23 

1 

2 

0.10 

5 

3 

0.10 

10 

4 

0.07 

15 

5 

0.07 

20 

6 

0.07 

25 

7 

0.07 

30 

8 

0.07 

35 

9 

0.065 

40 

10 

0.065 

45 

11 

0.06 

50 

12 

0.05 

55 

13 

0.05 

60 

Table  10  Source  function  C(xk)  calculated  according  to  solidification-boundary 


conditions  given 
14  to  17  (WELD 


in  Table  2,  where  A /  =  (12/60)  mm,  Xk=yk  =  0.0  for  k  =  1  to  13  and  z,k 
2). 
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15 


16 
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k 

cy  11.5 

Zk  (A  / ) 

1 

0.10 

1 

2 

0.10 

5 

3 

0.10 

10 

4 

0.10 

15 

5 

0.10 

20 

6 

0.10 

25 

7 

0.10 

30 

8 

0.10 

35 

9 

0.10 

40 

10 

0.10 

45 

11 

0.10 

50 

12 

0.06 

55 

13 

0.09 

60 

C, 


2.4 


2.4 


2.4 


2.4 


Xk  (  A/  ) 


-10.0 


10.0 


0.0 


0.0 


y*  (  A/) 


0.0 


0.0 


-10.0 


10.0 


constraint 
1  for  k  = 


10 


Table  11  Source  function  C(xk)  calculated  according  to  solidification-boundary  constraint 
conditions  given  in  Table  3,  where  Al  =  (15/60)  mm,  xu  =  yk=  0.0  for  k  =  1  to  13  (WELD  3). 


k 

cy  15. 5 

Zk  ( A/ ) 

1 

0.40 

1 

2 

0.10 

5 

3 

0.10 

10 

4 

0.10 

15 

5 

0.10 

20 

6 

0.10 

25 

7 

0.10 

30 

8 

0.10 

35 

9 

0.10 

40 

10 

0.10 

45 

11 

0.08 

50 

12 

0.07 

55 

13 

0.04 

60 

Table  12  Source  function  C(xk)  calculated  according  to  solidification-boundary  constraint 
conditions  given  in  Table  4,  where  Al  =  (15/60)  mm,  Xk  =  yk  =  0.0  for  k  =  1  to  13  (WELD  4). 


k 

cy  15. 5 

Zk  (A  / ) 

i  1 

0.40 

1 

2 

0.10 

5 

3 

0.10 

10 

4 

0.10 

15 

5 

0.10 

20  | 

6 

0.10 

25 

7 

0.10 

30 

8 

0.10 

35 

9 

0.10 

40 

10 

0.10 

45 

11 

0.08 

50 

12 

0.08 

55 

13 

0.05 

60 
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Table  13  Source  function  C(xk)  calculated  according  to  solidification-boundary  constraint 
conditions  given  in  Table  5,  where  A/  =  (15/60)  mm,  Xk  =  yk  =  0.0  for  k  =  1  to  13  (WELD  5). 


k 

cy  17. 5 

Zk  (A  l) 

1 

0.40 

1 
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0.30 

5 
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0.30 

10 

4 

0.10 

15 
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0.10 

20 
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0.10 
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0.10 
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0.10 
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0.10 

40 
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0.10 

45 
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0.09 

50 
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0.09 

55 

13 

0.05 

60 

Table  14  Source  function  C(xk)  calculated  according  to  solidification-boundary  constraint 
conditions  given  in  Table  6,  where  M  =  (3.0/60)  mm,  Xk  =  yk  =  0.0  for  k  =  1  to  13  (WELD  6). 
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Zk  (A  l) 
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45 

11 
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50 

12 

0.09 

55 

13 

0.05 

60 
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Table  15  Source  function  C(xk)  calculated  according  to  solidification-boundary  constraint 
conditions  given  in  Table  7,  where  A l  =  (18/60)  mm,  x/(  =  >’/,  =  0.0  for  k  =  1  to  13  (WELD  7). 


k 

cy  17. 5 

Zk  (A  / ) 
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1.00 
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0.30 
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3 

0.20 

10 

4 

0.10 

15 
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0.10 

20 

6 

0.15 

25 

7 

0.15 

30 
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0.20 

35 

9 

0.10 

40 
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45 

11 

0.12 

50 

12 

0.15 

55 

13 

0.05 

60 

Table  16  Source  function  C(xk)  calculated  according  to  solidification-boundary  constraint 
conditions  given  in  Table  8,  where  Al  =  (20/60)  mm,  Xk  =  yk  =  0.0  for  k  =  1  to  13  (WELD  8). 
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Zk  (A  l) 
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0.40 
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0.38 
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0.35 
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0.15 
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0.15 
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Fig.  3  Experimentally  measured  transverse  weld  cross  sections  of  solidification  and  HAZ-edge 
boundaries  for  steel  laser  weld  (Weld  1). 


>  ■  U  Hi  I  ivn 

Fig.  4  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  1  for  solidification  boundary  (Weld  1). 
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t  =  20  A? 


t  =  30  Af 


T(°C) 


Fig.  5  Temperature  history 
information  given  in  Table 
=  0.8  m/min  (Weld  1). 


§  § 


(°C)  of  transverse  cross  section  of  weld  calculated  using  cross  section 
1  for  solidification  boundary,  where  At  =  Al/\,  Al  =  (12/60)  mm  and  V 
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Fig.  6  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  at  estimated  HAZ-edge  boundary 
calculated  using  constraints  defined  by  Eq.(2),  and  given  in  Table  1  for  solidification  boundary 
(Weld  1). 


t  =  0  t  =  10  A? 


t  =  40  At  t  =  50Af 


t  =  80  At  t  =  90A  t 
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Fig.  7  Temperature  history  (°C)  at  estimated  HAZ-edge  boundary  of  transverse  cross  section  of  weld 
calculated  using  constraints  defined  by  Eq.(2),  and  given  in  Table  1  for  solidification  boundary, 
where  \t  =  A//V,  Al  =  (12/60)  mm  and  V  =  0.8  m/min  (Weld  1). 
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Fig.  8  Experimentally  measured  transverse  weld  cross  sections  of  solidification  and  HAZ-edge 
boundaries  for  steel  laser-GMA  hybrid  weld  (Weld  2). 


T(°  C) 

Fig.  9  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  2  for  solidification  boundary  (Weld  2). 
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t  =  40  At 
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Fig.  10  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross  section 
information  given  in  Table  2  for  solidification  boundary,  where  Af  =  A//V,  M  =  (12/60)  mm  and  V 
=  1.4  m/min  (Weld  2). 
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■  M  ^lr(°C) 

Fig.  11  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  at  estimated  HAZ-edge  boundary 
calculated  using  constraints  defined  by  Eq.(2),  and  given  in  Table  2  for  solidification  boundary 
(Weld  2). 
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Fig.  12  Temperature  history  (°C)  at  estimated  HAZ-edge  boundary  of  transverse  cross  section  of 
weld  calculated  using  constraints  defined  by  Eq.(2),  and  given  in  Table  2  for  solidification  boundary, 
where  \t  =  A//V,  A l  =  (12/60)  mm  and  V  =  1.4  m/min  (Weld  2). 
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Fig.  13  Experimentally  measured  transverse  weld  cross  sections  of  solidification  and  HAZ-edge 
boundaries  for  steel  laser-GMA  hybrid  weld  (Weld  3). 
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Fig.  14  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  3  for  solidification  boundary  (Weld  3). 
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Fig.  15  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross  section 
information  given  in  Table  3  for  solidification  boundary,  where  At  =  A//V,  Al  =  (15/60)  mm  and  V  = 
1.1  m/min  (Weld  3). 
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Fig.  16  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  at  estimated  HAZ-edge  boundary 
calculated  using  constraints  defined  by  Eq.(2),  and  given  in  Table  3  for  solidification  boundary 
(Weld  3). 
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Fig.  17  Temperature  history  (°C)  at  estimated  HAZ-edge  boundary  of  transverse  cross  section  of 
weld  calculated  using  constraints  defined  by  Eq.(2),  and  given  in  Table  3  for  solidification  boundary, 
where  At  =  A//V,  A /  =  (15/60)  mm  and  V  =  1.1  m/min  (Weld  3). 


Laser-GMA 
Hybrid  Weld 


Solidification  .  * 
Boundary  •  * 


Grade  S235 
Structural  Steel 


,  •  Laser  Power  15  kW 
.  •  Arc  Power  3.6  kW 

•  • 

.  Weld  Speed  1 .3  m/min 

] .  Laser-GMA  Gap  0 .5  mm 


Estimated  HAZ 


edge 


Fig.  18  Experimentally  measured  transverse  weld  cross  sections  of  solidification  and  HAZ-edge 
boundaries  for  steel  laser-GMA  hybrid  weld  (Weld  4). 
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Fig.  19  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  4  for  solidification  boundary  (Weld  4). 
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Fig.  20  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross  section 
information  given  in  Table  4  for  solidification  boundary,  where  A t  =  A l /V,  A/  =  (15/60)  mm  and  V  = 
1.3  m/min  (Weld  4). 
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Fig.  21  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  at  estimated  HAZ-edge  boundary 
calculated  using  constraints  defined  by  Eq.(2),  and  given  in  Table  4  for  solidification  boundary 
(Weld  4). 


31 


f-0 


t  =  20AJ 


t  =  40  Af 


Wl„ 

[AZ  Edge  1^', 

1 

u 

t  =  60  At 


t  =  30Ar 


t  =  70  At 


32 


t  =  80  At 


Fig.  22  Temperature  history  (°C)  at  estimated  HAZ-edge  boundary  of  transverse  cross  section  of 
weld  calculated  using  constraints  defined  by  Eq.(2)  given  in  Table  4  for  solidifcation  boundary, 
where  At  =  A//V,  A l  =  (15/60)  mm  and  V  =  1.3  m/min  (Weld  4). 
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Fig.  23  Experimentally  measured  transverse  weld  cross  sections  of  solidification  and  HAZ-edge 
boundaries  for  steel  laser-GMA  hybrid  weld  (Weld  5). 
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Fig.  24  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  5  for  solidification  boundary  (Weld  5). 
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Fig.  25  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross  section 
information  given  in  Table  5  for  solidification  boundary,  where  At  =  A//V,  A l  =  (15/60)  mm  and  V  = 
0.9  m/min  (Weld  5). 
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Fig.  26  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  at  estimated  HAZ-edge  boundary 
calculated  using  constraints  defined  by  Eq.(2),  and  given  in  Table  5  for  solidification  boundary 
(Weld  5). 
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Fig.  27  Temperature  history  (°C)  at  estimated  HAZ-edge  boundary  of  transverse  cross  section  of 
weld  calculated  using  constraints  defined  by  Eq.(2)  given  in  Table  5  for  solidifcation  boundary, 
where  At  =  A//V,  A l  =  (15/60)  mm  and  V  =  0.9  m/min  (Weld  5). 
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Fig.  28  Experimentally  measured  transverse  weld  cross  sections  of  solidification  and  HAZ-edge 
boundaries  for  steel  laser-GMA  hybrid  weld  (Weld  6). 
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Fig.  29  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  6  for  solidification  boundary  (Weld  6). 
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Fig.  30  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross  section 
information  given  in  Table  6  for  solidification  boundary,  where  At  =  A//V,  A l  =  (15/60)  mm  and  V  = 
0.9  m/min  (Weld  6). 
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Fig.  31  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  at  estimated  HAZ-edge  boundary 
calculated  using  constraints  defined  by  Eq.(2),  and  given  in  Table  6  for  solidification  boundary 
(Weld  6). 
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Fig.  32  Temperature  history  (°C)  at  estimated  HAZ-edge  boundary  of  transverse  cross  section  of 
weld  calculated  using  constraints  defined  by  Eq.(2)  given  in  Table  6  for  solidifcation  boundary, 
where  At  =  A//V,  A l  =  (15/60)  mm  and  V  =  0.9  m/min  (Weld  6). 
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Fig.  33  Experimentally  measured  transverse  weld  cross  sections  of  solidification  and  HAZ-edge 
boundaries  for  steel  laser-GMA  hybrid  weld  (Weld  7). 


Fig.  34  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  7  for  solidification  boundary  (Weld  7). 
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Fig.  35  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross  section 
information  given  in  Table  7  for  solidification  boundary,  where  At  =  A//V,  A /  =  (18/60)  mm  and  V  = 
0.9  m/min  (Weld  7). 
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Fig.  36  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  at  estimated  HAZ-edge  boundary 
calculated  using  constraints  defined  by  Eq.(2),  and  given  in  Table  7  for  solidification  boundary 
(Weld  7). 
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Fig.  37  Temperature  history  (°C)  at  estimated  HAZ-edge  boundary  of  transverse  cross  section  of 
weld  calculated  using  constraints  defined  by  Eq.(2)  given  in  Table  7  for  solidifcation  boundary, 
where  At  =  A//V,  A /  =  (18/60)  mm  and  V  =  0.9  ni/min  (Weld  7). 
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Fig.  38  Experimentally  measured  transverse  weld  cross  sections  of  solidification  and  HAZ-edge 
boundaries  for  steel  laser-GMA  hybrid  weld  (Weld  8). 
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Fig.  39  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  8  for  solidification  boundary  (Weld  8). 
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Fig.  40  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross  section 
information  given  in  Table  8  for  solidification  boundary,  where  At  =  A//V,  A l  =  (20/60)  mm  and  V  = 
0.5  m/min  (Weld  8). 
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Fig.  41  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  at  estimated  HAZ-edge  boundary 
calculated  using  constraints  defined  by  Eq.(2),  and  given  in  Table  8  for  solidification  boundary 
(Weld  8). 
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Fig.  42  Temperature  history  (°C)  at  estimated  HAZ-edge  boundary  of  transverse  cross  section  of 
weld  calculated  using  constraints  defined  by  Eq.(2)  given  in  Table  8  for  solidifcation  boundary, 
where  At  =  A//V,  A /  =  (20/60)  mm  and  V  =  0.5  m/min  (Weld  8). 

Discussion 

The  inverse  analysis  procedure  entails  calculating  a  three-dimensional  temperature  field  using 
experimentally  measured  constraint  conditions  on  the  solidification  boundary  This  temperature  field 
permits  calculation  of  temperature  histories  as  a  function  of  transverse  position  within  the  cross 
section  of  the  weld.  Shown  in  Figs.  4,  9,  14,  19,  24,  29,  34  and  39  are  two-dimensional  slices  of  the 
calculated  three-dimensional  temperature  fields  obtained  using  the  constraint  conditions  given  in 
Tables  1  through  8  for  the  measured  solidification  boundary,  which  are  parallel  to  the  relative 
motion  of  laser  or  laser-GMA  source  and  workpiece.  Shown  in  Figs.  5,  10,  15,  20,  25,  30,  35  and  40 
are  two-dimensional  slices  of  this  three-dimensional  temperature  field  that  are  perpendicular  to  the 
relative  motion  of  laser  or  laser-GMA  source  and  workpiece.  Referring  to  these  figures,  it  should  be 
noted  that  t  =  0  has  been  assigned  arbitrarily  to  a  two-dimensional  slice  at  the  leading  edge  of  the 
solidification  boundary.  Accordingly,  shown  in  these  figures  is  passage  with  time  of  the  calculated 
three-dimensional  solidification  boundaries  through  experimentally  measured  transverse  cross 
sections  of  these  boundaries,  which  are  indicated  by  sparse  dotted  contours.  For  the  planar  slices  of 
the  calculated  temperature  fields  shown  in  Figs.  4,  5,  9,  10,  14,  15,  19,  20,  24,  25,  29,  30,  34,  35,  39 
and  40,  the  weld  meltpool  and  regions  having  temperatures  below  melting  are  indicated  by  uniform 
black  and  banded  gray  scale,  respectively.  Referring  to  the  calculated  temperature  fields  shown  in 
these  figures,  the  constraint  conditions  on  the  calculated  three-dimensional  solidification  boundaries 
are  such  that  projections  of  all  their  two-dimenional  transverse  slices,  as  a  function  of  time,  are 
consistent  with  the  experimentally  measured  transverse  cross  sections  of  these  boundaries. 

Shown  in  Figs.  6,  11,  16,  21,  26,  31,  36  and  41  are  two-dimensional  slices  of  the  calculated 
three-dimensional  temperature  field  at  estimated  HAZ-edge  boundaries  obtained  using  the  constraint 
conditions  defined  by  Eq.(2),  and  given  in  Tables  1  through  8  for  the  solidification  boundary,  which 
are  parallel  to  the  relative  motion  of  laser  or  laser-GMA  source  and  workpiece.  Shown  in  Figs.  7,  12, 
17,  22,  27,  32,  37  and  42  are  two-dimensional  slices  of  this  three-dimensional  temperature  field  that 
are  perpendicular  to  the  relative  motion  of  laser  or  laser-GMA  source  and  workpiece.  Again, 
referring  to  these  figures,  it  should  be  noted  that  t  =  0  has  been  assigned  arbitrarily  to  a  two- 
dimensional  slice  at  the  leading  edges  of  the  isothermal  boundaries.  Accordingly,  shown  in  these 
figures  is  passage  with  time  of  the  calculated  three-dimensional  isothermal  boundaries  through 
experimentally  measured  transverse  cross  sections  of  the  estimated  HAZ-edge  boundaries,  which  are 
indicated  by  sparse  dotted  contours.  For  the  planar  slices  of  the  calculated  temperature  fields  shown 
in  Figs.  6,  7,  11,  12,  16,  17,  21,  22,  26,  27,  31,  32,  36,  37,  41  and  42,  regions  having  temperatures 
below  and  above  calculated  isothermal  boundaries  at  estimated  HAZ  edges  are  indicated  by  uniform 
black  and  banded  gray  scale,  respectively.  Referring  to  the  calculated  temperature  fields  shown  in 
these  figures,  it  can  be  seen  that  the  predicted  temperature  histories  at  cross  section  locations  close  to 
the  experimentally  estimated  HAZ  edge  are  within  a  range  of  temperatures  characteristic  of  that  for 
steels  (see  Fig.  1).  That  is  to  say,  these  results  show  reasonable  consistency  between  model  input. 
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solidification  cross  section  measurements,  and  model  output,  predicted  temperature  histories  close  to 
the  HAZ  edge. 

A  fundamental  aspect  of  inverse  analysis  methods  is  that  errors  associated  with  experimental 
measurements,  which  are  adopted  as  constraints,  are  encoded  onto  the  optimized  values  of  model 
parameters.  Accordingly,  quantitative  inverse  analyses  should  consider  the  sensitivity  of  calculated 
field  quantities  with  respect  to  measurement  errors  that  could  be  associated  with  constraint 
conditions.  It  follows  that  the  predicted  temperature  histories  at  the  HAZ  edges  are  expected  to  have 
errors  due  to  the  fact  that,  in  practice,  it  is  difficult  to  measure  solidification  and  HAZ  edge 
boundaries  for  steels.  In  particular,  metallographic  analysis  of  weld  cross  sections  [38],  using 
etchants,  must  consider  the  finite  thickness  of  the  mushy  zone  adjacent  to  the  solidus  isotherm,  as 
well  as  the  many  different  crystallographic  zones  in  proximity  of  the  HAZ  edge  (see  Fig.l),  which  in 
principle  may  not  be  well  defined.  Thus,  the  inverse  analysis  procedure  provides  a  consistency  check 
with  respect  to  experimental  procedures  for  measurements  of  solidification  and  HAZ  boundaries,  for 
a  given  steel  weld. 

The  results  of  this  study  can  be  adopted  for  more  efficient  inverse  thermal  analysis  of  other  types 
of  stainless  steel  welds.  This  follows  in  that  parameter  optimization  can  be  made  more  efficient 
using  initial  estimates  of  the  parameter  values,  which  require  only  fine  adjustment  with  respect  to 
constraint  conditions.  Model-parameter  optimization  for  welds,  whose  process  parameters  are  within 
similar  regimes  to  those  for  which  model-parameter  values  have  been  determined  previously  (e.g., 
parameters  C(xk),  xk ,  At  k  and  Tm  for  the  parametric  model  defined  by  Eqs.  (l)-(4)),  can  adopt 
these  values  as  initial  estimates  for  subsequent  fine  adjustment.  In  addition,  although  the  thermal 
diffusivity  k  and  melt  temperature  Tm  of  steels  may  vary,  this  variation  is  not  over  a  wide  range  of 
values.  This  is  the  case  in  general  for  different  types  metals  and  their  alloys.  It  follows  that 
parameter  optimization  for  a  specific  type  of  steel  weld,  which  uses  initial  estimates  of  parameter 
values  corresponding  to  different  types  of  steel  welds,  can  adopt  k  and  Tm  as  adjustable  parameters, 
as  well  as  other  model  parameters,  e.g.,  C(xk),  xk,  At.  Accordingly,  the  parametric  temperature 
histories  constructed  according  to  Tables  9-16  can  contribute  to  a  parameter  space  containing  a 
sufficient  range  of  parameters  corresponding  to  different  welding  processes,  process  conditions  and 
different  types  of  metals  and  their  alloys.  It  follows  that,  given  a  sufficient  accumulation  of 
parameterized  temperature  histories,  spanning  a  wide  range  of  process  conditions,  further 
investigation  should  concern  determination  of  an  optimal  structure  for  such  a  parameter  space. 

Finally,  which  should  be  acknowleged,  there  exists  many  different  types  of  steels,  for  whose 
welds  measurement  of  solidification  and  HAZ  edge  boundaries  may  be  difficult  and  possibly  not 
convenient.  A  classification  of  steels  that  can  be  characterized  for  quantitative  estimation  of 
solidification  and  HAZ  edge  boundaries  is  beyond  the  scope  of  this  study.  The  present  study 
indicates,  however,  that  for  steels  permitting  quantitative  estimation  of  solidification  and  HAZ  edge 
boundaries,  temperature  histories  can  be  calculated,  which  are  consistent  with  bounding 
temperatures  of  the  HAZ.  This  consistency  provides  a  validation  of  model  parameter  values. 

Conclusion 

A  specific  objective  of  this  report  is  to  further  examine,  for  the  case  of  structural  steel  welds,  the 
concept  of  using  experimentally  measured  temperature-field  constraints  for  inverse  thermal  analysis. 
In  particular,  examination  of  a  procedure  for  accessing  the  consistency  of  calculated  temperature 
histories  for  steel  welds  that  uses  experimentally  measured  estimates  of  the  HAZ  edge.  A  general 
objective  of  this  report  is  calculation  of  parametric  temeprature  histories  by  inverse  thermal  analysis 
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for  structural  steel  deep-penetration  welds  corresponding  to  various  weld  process  parameters.  This 
report  contributes  to  the  continuing  evolution  of  a  parametric  representation  of  the  temperature  field 
for  inverse  thermal  analysis  of  welds  associated  with  different  types  of  metals,  their  alloys  and  weld 
process  conditions 
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